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Monte Carlo simulations of slow positron production via energetic electron interaction with a solid target 
have been performed. The aim of the simulations was to determine the expected slow positron beam intensity 
from a low energy, high current electron accelerator. By simulating (a) the fast positron production from a 
tantalum electron-positron converter and (b) the positron depth deposition profile in a tungsten moderator, 
the slow positron production probability per incident electron was estimated. Normalizing the calculated 
result to the measured slow positron yield at the present AIST LINAC the expected slow positron yield as a 
function of energy was determined. For an electron beam energy of 5 MeV (10 MeV) and current 240 /iA (30 
/jA) production of a slow positron beam of intensity 5 x 10 6 s _1 is predicted. The simulation also calculates 
the average energy deposited in the converter per electron, allowing an estimate of the beam heating at a 
given electron energy and current. For low energy, high-current operation the maximum obtainable positron 
beam intensity will be limited by this beam heating. 



I. INTRODUCTION 

The use of electron accelerators is a well established 
technique for the production of slow positrons, with nu- 
merous facilities based on linear accelerators (LINAC) 1 - - — 
or microtron o 10 ' 11 . Positrons are generated via pair cre- 
ation when energetic electrons are stopped in a high-Z 
target, the electron-positron converter. These positrons 
are then moderated, i.e. slowed to thermal energies, via 
interaction with a suitable material such as tungsten, and 
a slow positron beam produced. The positron produc- 
tion probability increases from threshold (1 MeV) with 
increasing electron energy and typically a value between 
14 MeV and 100 MeV has been used. 

However, with the reduction in cost of low energy (< 15 
MeV), high-current industrial LINACs and development 
of new types of accelerator such as the Rhodotron 1 - 2 - it 
has become appealing to consider producing intense, slow 
positron beams using low energy, high current electron 
beams. One group, at Saclay, is developing such a system 
based on a 6 MeV commercial LINA C 13 ' 14 . 

For the past 20 years the AIST LINAC has been used 
to produce high intensity, slow positron beams for ma- 
terials research 4,15 . A schematic of the AIST electron- 
positron converter and moderator assembly is shown in 
ngureQ] Electrons with a kinetic energy of 70 MeV are di- 
rected onto a water cooled tantalum disk with a thickness 
of 6 mm. Positrons emerging in the forward direction 
are incident on a moderator composed of strips of 50 /im 
tungsten foil arranged in a rectangular mesh. Slow (mod- 
erated) positrons are electrostatically extracted from the 
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FIG. 1. A schematic diagram of the converter and moderator 
assembly currently in use at the AIST LINAC. Electrons at 70 
MeV are directed on to a 6 mm thick, water cooled tantalum 
block in air. Positrons emerging in the forward direction pass 
through a thin titanium film into vacuum and are incident on 
an array of 50 fim. tungsten films arranged in a rectangular 
mesh. Moderated positrons at low energy (several eV) are ex- 
tracted from the moderator assembly by the applied positive 
potential («10 V) and magnetically guided into a transport 
beamline. 



moderator then formed into a beam and magnetically 
guided to experimental stations located more than 20 m 
from (and well shielded from) the radiation produced in 
the converter. 

The AIST LINAC operates at a frequency of 100 
Hz, delivering a pulse of around 30 nC with a pulse 
width of 1/is. Slow positrons are produced with a 
similar time structure which is impractical for use in 
positron annihilation spectroscopy (PAS) type experi- 
ments such as Doppler broadening of annihilation ra- 
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diation and positron annihilation lifetime spectroscopy 
(PALS) 16 . In practice, the beam is stored in a linear 
trap and transformed into a quasi-DC beam. In PALS 
this DC beam is chopped and bunched to form a pulse 
train with a pulse width of around 100 ps and a frequency 
of around 10 MHz. 

This manipulation of the beam timing structure leads 
to a considerable loss of positrons, typically at most 10% 
reach the sample in a PALS experiment. If the electron 
accelerator delivered smaller electron pulses with much 
higher frequency this situation could be improved. This 
is possible with a superconducting accelerator (SCA) and 
our group is now developing a dedicated SCA for positron 
productio n 17 ' 18 . In an ideal case the positron beam pro- 
duced in the converter/moderator could be used in PALS 
experiments without manipulation, resulting in an in- 
creased beam transport efficiency. 

However, in comparison to the beam energy used at 
the current AIST LIN AC (70 MeV), the proposed SCA 
will be based on an accelerating module with a maximum 
acceleration of around 7.5 MV. Two such modules, pre- 
viously used on the JAERI FEL projec t 19 ' 20 , have been 
obtained by our group. Using these modules in a single 
pass configuration we expect a maximum beam energy 
of 7.5 MeV (1 module) or 15 MeV (2 modules). How- 
ever, initially we expect to run with slightly lower than 
maximal field gradient with an expected beam energy of 
around 5 MeV per module. 

Using an accelerator with a low beam energy sub- 
stantially reduces the positron production probability al- 
though it also leads to a much reduced radiation dose and 
induced activity from the accelerator. The maximum of 
the cross section for neutron production via the giant res- 
onance occurs around 14 MeV— for tantalum. Keeping 
the electron energy below this maximum reduces the re- 
quired radiation shielding dramatically, an advantage for 
compact, low-cost systems. 

The purpose of the present paper is to estimate the 
yield of slow positrons from the proposed SCA. We report 
results of Monte Carlo simulations of the converter and 
moderator using the Penelope200822 code. Penelope2008 
can perform calculations over a wide energy range, from 
a few hundred eV to about 1 GeV and provides full sup- 
port for positron interactions. It is thus highly suited to 
modeling of positron production from energetic electron 
beams. 

Previously, several groups have reported Monte 
Carlo simulations of positron production from electron 
accelerators 2 ^—. In particular Gagliardi and Hunt 28 
reported a systematic study of tungsten converters at 
normal and glancing irradiation for a range of incident 
electron energies. Most previous simulations have con- 
centrated on the fast positron production from the con- 
verter. In recent years there has been progress in devel- 
oping low energy extensions to the Monte Carlo code in 
order to fully simulate the moderation proces o 29 ' 30 . In 
the present study both the converter and moderator are 
simulated using a standard code (without low energy ex- 



tensions) by implementing a well established model of 
the moderation process. 

II. SIMULATIONS 

A diagram of the simulation geometry is shown in fig- 
ure [U The simulation is divided into two parts which 
are simulated separately, (a) the converter and (b) the 
moderator. For each simulation between 10 6 and 10 8 
primary particle trajectories are calculated, with larger 
simulations typically necessary at lower electron energies 
when positron production probabilities are reduced. The 
low energy cut-offs were typically set at 50 keV for elec- 
trons and positrons and 5 keV for photons. Reducing 
the cut-offs below these values caused the computation 
time to increase without any significant change in the 
simulation results. 



A. The electron-positron convertor 

In this simulation an electron beam is incident on a 
cylindrical tantalum target. Input parameters are the 
electron energy, E e , incident angle, a, and convertor 
thickness, I. The output parameters of most interest are 
the angular, j(8,<p), and energy, j(E p ) distributions of 
emerging positrons. Here E p is the energy of the emitted 
positron and 8, 4> the polar and azimuthal angles respec- 
tively. The polar angle is defined with respect to the 
axis perpendicular to the converter faces. The azimuthal 
angle is cf> = in the plane defined by the axis perpen- 
dicular to the converter face and the electron beam di- 
rection. When electrons are normally incident (a = 0) 
the simulation is symmetric in the angle <f>. Finally we 
also note the average energy deposited in the converter 
per electron, E$. 

1. Optimum converter thickness 

Initially the optimum converter thickness, l op t, was de- 
termined. At each set of input parameters the probability 
(per electron) of producing a fast positron which leaves 
the converter in the forward direction (0° < 6 < 90°) was 
calculated, 7/, and is given by, 

7/ = / / 2 7(6>, </>)sm0d0d</> (1) 
Jo Jo 

Figure [3] shows the results of a these simulations for 
a range of electron energies between E e = 5 MeV and 
100 MeV and converter thickness up to 10 mm. At each 
electron energy the positron production probability, 7/, 
increases to a maximum before decreasing with increas- 
ing converter thickness. This optimum thickness, l op t, in- 
creases with increasing E e and the peak becomes broader. 
Figure H] shows a plot of the optimal thickness, l opt as a 
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a) Converter 




b) Moderator 



FIG. 2. Schematic of the (a) converter and (b) moderator 
simulation geometry with calculated output parameters in 
boldface. In (a) electrons of energy E e , at an angle a to the 
normal, are incident on a cylindrical tantalum target with 
variable thickness, I, and the energy j(E p ), angular j(6,(J)) 
distributions of emerging positrons calculated. The average 
energy deposited in the converter per incident electron, Ed, is 
also calculated. In (b) mono-energetic positrons at an angle 
a p are incident on a 50 /im tungsten foil and the positron 
depth deposition distribution, f(z), determined. Using this 
distribution the probability of re-emission of a slow positron, 
j), is calculated as outlined in the text. By convolving the 
results of simulation (a) and (b) , the probability of producing 
a slow (moderated) positron per incident electron, j a , can be 
estimated. 



function of electron energy, E e , for electrons at both nor- 
mal (a — 0°) and glancing (a = 87°) incidence. The 
error bars represent not the statistical error in the cal- 
culation but rather the range over which 7/ is > 90% of 
the maximum. 



10 ■ 



10' 



10 ; 



10" 



AAA A A 



-*■ — * — ► 100 

♦ ♦ 40 



-* 20 
* 15 

-• 10 



5MeV 



4 6 8 

Converter Thickness, / (mm) 



10 



FIG. 3. (Color online) Variation in positron production prob- 
ability in the forward direction, 7/, with converter thickness, 
I, for a range of electron energies between E e = 5 MeV and 
100 MeV. Electrons are normally incident on a tantalum con- 
verter in every case. The dashed line connects values at the 
optimal converter thickness, l op t, at each electron energy. 
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FIG. 4. (Color online) Optimum (maximum probability of 
fast positron production in the forward direction, 7™ aa: ) con- 
verter thickness, l opt , calculated as a function of electron en- 
ergy, Ee., at both normal (a = 0°) and glancing (a = 87°) 
incidence. The error bars represent the range of thickness 
over which 7/ > 0.9 r yJ lax . Also shown is an extrapolation 
based on previous measurements by Ley^i and the model of 
Akahane et aZ.— . 



At normal incidence, for the range of energies studied of electron energyj E ^ as; 



lopt varies from 0.6 mm at 5 MeV to 7.5 mm at 100 
MeV. The calculated value at 70 MeV is 6.5 mm, close 
to the converter currently in use at AIST (6 mm). For 
glancing irradiation there is almost no variation in l opt 
with E e , l opt changes from 0.4 mm at 5 MeV to 0.85 mm 
at 100 MeV. In all cases the value of l op t quoted is the 
real thickness of the converter material, i.e. not scaled 
by the angle of irradiation, a. For normal incidence, a fit 
to these calculated data points allows us to estimate the 
optimum thickness, l op t, of a Ta converter as a function 



l opt [mm] = 0.6 + 0.64(£: e [MeV] - 5) 



0.54 



(2) 



Also shown on figure [4] is an extrapolation based on 
measurements of the optimum converter thickness at var- 
ious laboratories 31 , ! opi [iran] = 0.67 + 0.0953(S e [MeV]). 
This linear extrapolation tends to over-estimate the op- 
timal thickness at high energy. Better agreement is 
seen with the simple model developed by Akahane et 
al. which assumes that the initial energy of the electron 
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is shared between electron-positron and bremsstrahlung 
production^. They estimate an optimal thickness for any 
material based on the radiation length, XQ, as; 

l opt = X0^, (3) 

ti c 

where E c is the critical energy below which ionization be- 
comes the dominant electron energy-loss mechanism. For 
tantalum E e = 11 MeV and XO = 4.1 mm. This simple 
model provides a good estimate of l op t at high electron en- 
ergy but underestimates the optimum value for E e < 20 
MeV. However, for intermediate energies both the linear 
extrapolation and the Akahane model provide reasonable 
estimates of the optimal converter thickness. 

2. Variation of fast positron production with electron 
energy 

Assuming optimal thickness converters are used we can 
then plot the positron production probability in the for- 
ward direction, 7/, as a function of electron energy, E e . 
Figure [5] shows this plot for normally incident electrons 
(a = 0°) and shows a reduction in 7/ by 3 orders of mag- 
nitude when E e is reduced from 70 to 5 MeV. Analysis 
for electrons incident at a = 45° and 87° is also shown 
on figure [5] The values are plotted as a comparison to 
the normally incident beam, i.e. 7/[ct]/7/[a = 0°]. In all 
cases the maximum value of 7/ using converters of opti- 
mal thickness are plotted. It is clear that while the opti- 
mal converter thickness l op t is sensitive to the angle of ir- 
radiation (figure 2]), the maximum available fast positron 
intensity is less sensitive, provided optimized converters 
are used. At all energies 7/ is increased when a = 45° 
with a maximum increase of 13% at for E e = 5 MeV. 
Glancing irradiation (a — 87°) reduces 7/ at low elec- 
tron energies but actually leads to increased production 
for energies > 70 MeV. 

3. Positron energy and angular distributions 

The calculated energy distribution of positrons emerg- 
ing in the forward direction, r yf(E p ) 1 is plotted in figure[S] 
for E e = 5 MeV and 70 MeV for a = 0°, 45° and 87°. The 
energy distributions are very broad, covering the whole 
energy range from the low energy cut-off at 50 keV up to 
almost the initial electron energy E e . For E e = 70 MeV 
the distribution is peaked around 2-3 MeV while at E e 
= 5 MeV the peak is around 0.6-0.7 MeV. 

Figure [7] plots the polar angular distribution of emerg- 
ing positrons in the tf> = plane after irradiation with (a) 
70 MeV and (b) 5 MeV electrons at a = 0°, 45° and 87°. 
For non-normal irradiation the emerging positron distri- 
bution depends on both the polar (9) and azimuthal (</>) 
angles. At E e — 70 MeV, the polar distribution is peaked 
in the same direction as the incident electron beam for a 
— 0° and a = 45°. For a — 87° the peak is broader and 



lO" 1 ^ 




A 

A 

* A 


1 1 1 1 1 1 1 


* 

A 


10 ' 2 " 








* 




10-*- 




/ * 

/ * 
*/ * 


— "— Y ( [a = 0°] 
* Y( [a = 45°]/7 f [a = 0°] 










* 7, [a = 87°] / Y f [a = 0°] 




( 


) 


20 


1 1 1 1 1 1 1 

40 60 80 
Electron Energy, E (MeV) 


100 



FIG. 5. (Color online) Fast positron production probability 
in the forward direction, 7/ at I = l op t, calculated for sev- 
eral input electron energies E e for electrons normally incident 
(a = 0°) incident on the converter. Also shown in the rela- 
tive increase or decrease in 7/ for a — 45° and 87°. Tantalum 
converters of optimal thickness were used for each simulation. 




Positron Energy, E (MeV) 

FIG. 6. (Color online) Typical energy distributions of 
positrons emerging in the forward direction, 7/(-B p ), after ir- 
radiation with 5 MeV and 70 MeV electrons at a = 0°, 45° 
and 87°. In all cases tantalum converters of optimum thick- 
ness, lopt, were simulated. 



has a maximum around 9 = 70° and there is significant 
positron production in the backward direction, actually 
greater than the forward production in this case. 

For low energy irradiation (E e — 5 MeV) the polar dis- 
tribution is much broader and peaked around the normal 
(9 = 0) for all incident angles. Although not plotted, the 
azimuthal distribution is also more uniform at E e = 5 
MeV compared to the case for E e = 70 MeV where the 
emerging positrons are emitted in a smaller cone around 
the axis. 

For normally incident electrons, a = 0, the emerging 
positron angular distribution 7(0, </)) is symmetric about 
the azimuthal angle cj) so 7(6*, <j>) can be described by the 
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FIG. 8. (Color online) Angular distribution of positrons in the 
forward direction from a Ta converter of optimal thickness 
when irradiated at normal incidence by electrons of energy 
E e . The angular distributions were fitted with the function 
7(0) = no[exp(9/9 c )—A] and values of characteristic angle, 9 C , 
and constant, A, are plotted. As the electron energy increases 
the positrons tend to emerge closer to the axis (9 = 0). Also 
shown is the angle 6^ where the intensity falls to half that at 
9 = 0, i.e. j(0 h ) = 0.5 n . 



l(Qh) = 0.5 no. This value is also plotted on figure[5]and 
shows a decrease from 46° at 5 MeV to 13° at 100 MeV. 
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FIG. 7. (Color online) Typical polar angular distributions 
in the central plane, j(9,cj) = 0), of emerging positrons after 
irradiation with (a) 70 MeV and (b) 5 MeV electrons at a — 
0°, 45° and 87°. In all cases tantalum converters of optimum 
thickness, l op t, were simulated. 



polar angle alone, 7(6*). For all energies the distribu- 
tions are peaked on axis (6 = 0) and have minima at 
9 = 90°. For 8 > 90° the positrons are emerging from 
the back of the converter in the opposite direction to 
the electron beam. This fraction is small at high en- 
ergy but is comparable to the production in the forward 
direction for E e = 5 MeV. For high energy irradiation 
it was calculated previously that high energy positrons 
should be emitted with an angular distribution of the 
forrn^ 70 = rig exp(6/8 c ), where no is the intensity on 
axis (9 = 0), and 9 C is a characteristic decay angle. This 
function was found to be a good description of 7(0) at 100 
MeV, but, for lower energies it was found to be necessary 
to introduce an offset in order to maintain an accurate 
fit, i.e. 7(0) = n [exp(9/9 c ) - A]. 

Figure [5] shows the fitting parameters for the angular 
distribution as a function of E e . As E e increases the 
offset tends to zero and 9 C decreases to around 20°. A 
second measure of the angular spread can be defined as 
the angle where the intensity drops to half that on axis, 



The previous Monte Carlo simulations can help us es- 
timate the total number and energy, angular distribu- 
tions of positrons leaving the converter. In practice these 
positrons are then moderated, a process where the ener- 
getic positron is implanted into a material, thermalizes 
and may, with some probability diffuse to the surface and 
be re-emitted with a low (several eV) energy. At present 
we use an array of 50 /im tungsten films arranged in a 
rectangular mesh (figure [I]) . 

Since it is impossible to simulate the full thermaliza- 
tion and diffusion process with the Monte Carlo code, a 
simple model is used as a basis for the present calcula- 
tion of moderator efficiency. Typically the moderation 
process is considered on the basis of a diffusion length, 
L+, and re-emission branching ratio, yo, for the mod- 
erator material. Since the probability that a positron 
can diffuse back to the surface from a given depth, z, 
is equal to exp(— z/L+y&, then the total probability of 
re-emission, 77, is given by; 

f-OC 

V = Vo f(z)cxp(-z/L + )dz (4) 
Jo 

where f(z) is the positron deposition depth profile. 

For polycrystalline tungsten, Suzuki et al. 3i measured 
a diffusion length, L + = 55 nm, and re-emission branch- 
ing ratio, yo = 0.27. A schematic of the simulation ge- 
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FIG. 9. (Color online) Depth penetration profile, f(z), of 
0.5 MeV positrons incident at a p =45° on a 50 fim W foil. 
The solid line shows f(z) scaled by the function exp (—z/L+) 
which is equal to the probability of a positron diffusing back 
to the surface. The total area under this curve scaled by the 
re-emission branching ratio, yo, determines the re-emission 
probability, 77. The inset shows the full depth penetration 
profile, f(z). 

ometry is shown in figure [2jb) . The implantation depth 
profile, f(z), of positrons incident on W at a given en- 
ergy, E p , and angle, a p , is calculated by the simulation 
and values of L + and yo quoted above used to determine 
the re-emission probability, 77, according to equation @] 
Re-emission from both the front and back faces of the 
50 /jm foil was included in the analysis. An example of 
the calculated depth profile, f(z), and the probability 
of diffusion to the W surface is shown in figure |H1 The 
re-emission probability is determined by calculating the 
total area under the solid curve, scaled by the re-emission 
branching ratio. 

Results of this analysis are shown in figure QJJ with the 
re-emission probability, rj, plotted as a function of inci- 
dent positron energy, E p , for a range of incident angles, 
a p . It is clear that positrons incident on the modera- 
tor with low energy and at shallow angles have a much 
high re-emission probability than fast, normally incident 
particles. 

Taking the angular and energy distributions of fast 
positrons from the previous converter simulation we can 
then estimate the slow (moderated) positron production 
probability, 7 S , at each initial electron energy. We con- 
sider only the case where electrons are normally incident 
on the converter. As a first step all positrons emitted at 
a polar angle 8 > 30° from the normal are ignored, an as- 
sumption based on the geometry of the current converter 
and moderator. The ratio, R(30°), is defined as the ratio 
of fast positrons emitted with angles less than or equal to 
30° to the total forward (0-90°) production, i.e. R(30°) = 

Io /6 '~ff( e ) de / Io /2 ! 7/(0)d0- R(30°) decreases from 0.81 
at 100 MeV to 0.57 at 5 MeV, for low energy operation 
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FIG. 10. (Color online) Probability of positron re-emission, rj, 
from a thin W foil (50 (im) as a function of incident positron 
energy, E p , and angle of irradiation, a p . 

it is therefore more important to locate the moderator as 
close as possible to the converter. 

We also assume that all fast positrons from the con- 
verter with 8 < 30° are incident on a W moderator foil 
at an incident angle a p = 90 - 6 (since the converter 
and moderator faces are perpendicular, see figure [lj. Fi- 
nally it is assumed that all re-emitted positrons can be 
extracted to the slow positron beam. 

For the fast positrons with 8 < 30°, the fraction in each 
of the angular ranges, 8 = 0-10°, 10°-20° and 20°-30° was 
calculated and the energy distribution of this component 
(assumed to independent of angle) was convoluted with 
the moderation probability function at an angle of 85°, 
75° and 65° respectively and the contributions summed. 

The results of this conversion are plotted in figure 1111 
which shows the slow (moderated) positron production 
probability, j s , as a function of E e . Although the fast 
positron production probability is lower for lower energy 
incident electrons, the positrons created have a lower en- 
ergy distribution (figure^ and are hence moderated with 
higher efficiency than those created with a high energy 
beam. This result is somewhat negated by the broader 
angular distribution of the emitted fast positrons (figure 
[7]) , however on balance the efficiency of moderation tends 
to increase for lower energy irradiation. The moderator 
efficiency e at a given energy is defined as; 

and is also plotted on figure [IT] e decreases from 1.72 x 
10~ 4 at 5 MeV to 2.42 x 10~ 5 at 100 MeV. 

The slow positron production probability, 7^ , at 5 MeV 
is around 0.4% that at 70 MeV, an improvement of a fac- 
tor of 4.6 compared to the same ratio for un-moderated, 
fast positrons due to the increased moderation efficiency, 
e. The calculated fast and slow positron production prob- 
abilities along with other relevant data are summarized 
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FIG. 11. (Color online) Probability of slow positron produc- 
tion per incident electron, 7 S , as a function of electron energy, 
E e . The data points are based on the fast positron production 
probability, 7/, calculated in the first simulation multiplied 
by a conversion factor which accounts for the moderation ef- 
ficiency. Also shown is a fit to the calculations using a power 
law. 



TABLE I. Predicted fast, 7/, and slow, 7 a , positron produc- 
tion probabilities from an optimal thickness tantalum con- 
verter and normal incidence electron irradiation. The ratio 
of fast positrons produced within 30° of the normal is also 
shown, R(30°), along with the calculated moderation effi- 
ciency, e, based on the model described in the text. Numbers 
in round brackets indicate error estimates and square brackets 
order of magnitude. 
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1.09(05) [-7] 


9.48(31) 


-5] 


15 


2.3 


6.48(08) 


-3] 


0.61 


2.71(10) [-7] 


6.80(08) 


-5] 


20 


3.0 


1.38(01) 


-2] 


0.64 


4.88(12) [-7] 


5.39(04) 


-5] 


40 


5.0 


6.12(07) 


-2] 


0.71 


1.53(02) [-6] 


3.54(04) 


-5] 


70 


6.5 


1.56(01) 


-1] 


0.78 


3.27(07) [-6] 


2.71(02) 


-5] 


100 


7.5 


2.56(01) 


-1] 


0.81 


5.02(08) [-6] 


2.42(01) 


-5] 
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III. COMPARISON TO MEASURED SLOW POSITRON 
YIELDS 

The maximum slow positron intensity from the present 
AIST facility was measured previously to be around 5 x 
10 7 s" 1 at a beam energy of 70 MeV and current of 10 
/j,A. This corresponds to a production probability of 8.35 
x 10 -7 e + /e~. The simulation result at this energy is 3.27 
x 10~ 6 , or around 4 times the measured value. There 
are several factors which account for this difference in- 
cluding; 1. Our assumption that all positrons emitted 
at a polar angle less than 30° are incident on the mod- 
erator foil. Some fraction at very low angles may pass 
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FIG. 12. Simulation result normalized to the AIST measure- 
ment at 70 MeV plotted alongside previous measurements at 
LLNIj2 and Argonne^ and Mitsubishi Electric Corporation 5 . 
The LLNL and Argonne measurements tend to underesti- 
mate the positron production at low energy compared to the 
present simulation however much better agreement is seen 
with the 18 MeV measurement at Mitsubishi. 



through without meeting any of the W foils; 2. We also 
assume that every positron which is moderated can be 
extracted to the slow positron beamline. Clearly some 
fraction will not be guided out of the moderator assem- 
bly due to re-collisions with the W surface; 3. Degra- 
dation of the moderator due to electron induced defects 
and surface contamination. 

In order to compare the current simulation results to 
the measured value at 70 MeV the simulation result was 
normalized by scaling by a factor of 25.6%. This normal- 
ized result, 7tv, is plotted in figure [T2l along with several 
previously reported measurements 3 ^' ^ 



The normalized simulation result shows good agree- 
ment with the measured results at LLNL especially at 
high energy. The LLNL LINAC used a Ta convertor 
and W moderator assembly similar to that considered 
here. The measurements had a minimum energy of 
17.5 MeV and were taken with a fixed assembly, i.e. 
fixed convertor thickness of 5 mm, rather than an op- 
timized convertor. The authors assumed a simple linear 
electron energy relationship with a cut-off at 15 MeV, 
7 s [e+/10 8 e~] = 1.4(E e [MeV] - 15), (fit fromS 5 .). 

The current simulation suggests a much higher pro- 
duction probability at lower energy (E e < 20 MeV) 
than that predicted by the linear LLNL fit. This dis- 
agreement is also apparent in the low energy Argonne 
measurement :] 8 ' 36 . Much closer agreement is seen with 
the measurement at the Mitsubishi Electric Corporation 
LINAC 5 where a 5 mm Ta converter and W moderator 
were used and an estimated slow positron production ef- 
ficiency of 1.35 xl0~ 7 s _1 at 18 MeV was reported. 
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IV. BEAM HEATING OF THE CONVERTER 

Based on the previous calculations it is then straight- 
forward to estimate the expected slow positron intensity 
at a given electron beam energy and current. It is also 
clear that if intense slow positron beams are to be gen- 
erated using a low energy accelerator then a high beam 
current is required. However, at this stage we need to 
consider the energy deposited in the converter. The ra- 
tio of average energy deposited in the converter per elec- 
tron, Ed, divided by the electron energy E e is plotted in 
figure fT3Ta) as a function of the converter thickness, I, 
for a normally incident electron beam. This ratio deter- 
mines the fraction of the beam power deposited in the 
converter and increases with increasing converter thick- 
ness up to a maximum. For the range of energies studied 
this maximum is around 0.8, showing that even with very 
thick converters around 1 /5 of the input power is not de- 
posited in the converter but instead escapes in the form of 
backscattered and transmitted electrons, positrons, and 
high energy photons. 

The energy deposition in the converter can also be de- 
scribed by the specific energy loss parameter, AE, which 
is defined as, 



AE= ^[MeV/egcm -2 ] 
I p 



(6) 



where I is the converter thickness and p is the material 
density. AE is plotted in figure fT3l b). At low electron 
energy AE quickly reaches a maximum and subsequently 
decreases rapidly with increasing thickness. As the elec- 
tron energy is increased the thickness at which this max- 
imum is reached increases and is found to be similar the 
thickness at which positron production in the forward 
direction reaches a maximum, l m ax- At low energy the 
maximum value of AE is greater than 2 MeV/e - g cm~ 2 , 
decreasing to around 2 MeV/e~ g cm" 2 at intermediate 
energy, before increasing again for E e > 70 MeV. 

As the electron beam current increases the energy de- 
posited will eventually be sufficient to melt or vaporize 
the converter material. An estimate for the cooling power 
of a water cooled converter, P w , is given by the following 
equatior 



.37 



P 
- 1 a 



{T b -T w )2ir\l 
(0.5 + ln(r 2 /ri)) 



(7) 



where r 2 and r\ are the radius of the converter and elec- 
tron beam respectively, A is the thermal conductivity and 
I the converter thickness. Xj, is the maximum tempera- 
ture of the material inside the electron beam (on axis), 
while T w is the temperature of the copper block in con- 
tact with the water coolant. For tantalum p = 16.7 g/cm 3 
and A = 57.6 W/(m K) at 300 K. The thermal conduc- 
tivity of tantalum has been measured to increase with 
increasing temperature, up to a value of 61.4 W/(m K) 
at 1800 K2£ but this small temperature dependence was 
ignored in the subsequent analysis and the room temper- 
ature value used. 
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FIG. 13. (Color online) Plots of (a) the fractional energy ab- 
sorption in the converter, Ed/E e , and (b) the specific energy 
loss in the converter per incident electron, AE, for an electron 
beam at normal incidence impinging on a tantalum converter 
as a function of thickness, The dashed lines connect val- 
ues at the optimal converter thickness, l op t, at each electron 
energy. 



The converter material will also radiate heat, with the 
power dissipation, P r , given by, 

P r = 2ecr( f ' (T(r) 4 - T 4 )27rrdr) (8) 
Jo 

where To is the temperature of the surrounding envi- 
ronment and a is the Stcfan-Boltzmann constant (5.67 
xlO -8 J s" 1 m~ 2 K~ 4 ). e is the spectral emissivity 
which depends on the condition of the materials surface, 
for un-polished tantalum it was measured to be around 
0.32 over the temperature range 1100 - 2300 XM. T{r) 
is the radial temperature profile and is given by, 

T (r)=T b -£ 1 C {0<r< ri } 



4ttA/ rf 
P 



2n\l 



[0.5 + ln(r/r 2 )] {n < r < r 2 } (9) 



At equilibrium the power deposited in the converter 
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by the electron beam is equal to the sum of the power 
I 



The maximum electron 



where T m is the melting point of tantalum (3270 K) and 
C and D are constant equals to ( s+ETr^Tnll an< ^ ^ ea 
respectively. 

It is clear that if the electron accelerator can pro- 
vide sufficient electron current the maximum number of 
positrons at a particular electron energy, Y max , is equal 
to the product of the slow positron production probabil- 
ity, 7 S , and the maximum possible electron current, I ma x, 
i.e., 

-t^ T 7S j 7S /I Q\ 

I max — 7s Imax OC OC . . \ *-^) 

In practice the yield will be lower than Y max as it is 
unpractical to operate at I m ax- The effective maximum 
current will be some fraction of I ma x- 



V. EXPECTED BEAM INTENSITY 
A. Optimal thickness converters, l opt 

Figure [Til shows a 2D surface plot of the expected slow 
positron yield, Y, as a function of electron beam energy 
and current according to, Y = I ■ At each set of 
input parameters the temperature of the converter mate- 
rial inside the electron beam Tj, was also estimated using 
equation [TJ with T w = 350 K, T = 300 K, r x = 1 cm 
and T2 = 2 cm. The thickness of the converter was also 
varied with energy according to equation [2] 

A comparison of the present AIST LIN AC and possible 
future operation with a higher current lower energy SCA 
is shown in table [TTJ The present accelerator operates 
at 70 MeV and a maximum current of 10 fiA. The table 
also shows the results for a high current beam (200 /iA). 
With a 5 MeV electron beam a positron yield of 4.2 x 
10 6 s~ x is expected. A 200 fiA (beam power = I kW) 
beam is well within the capabilities of the proposed SCA, 
therefore the main limitation to increasing the positron 
yield is the power deposited in the Ta convertor. The 
simulations suggest that at these conditions the power 
deposited is around 488 W. Using equation[T0]for a water 
cooled converter with the dimensions given above (rj = I 
cm and r-i = 2 cm) the estimated maximum temperature 
the converter material inside the electron beam, T u , is 



dissipated by the water cooling and radiative losses. 

(10) 



(11) 
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FIG. 14. Surface plot showing the calculated slow positron 
yield, (Y), as a function of incident electron energy E e (x-axis) 
and current I e (y-axis). Ta converters of optimal thickness 
lopt and the normalized slow positron production probability 
7jv (figure I12p are used to calculate Y. The temperature of 
the converter material inside the electron beam Tb was also 
calculated according to equation [TU] and is also shown. The 
high-current cut-off is determined as the current at which the 
Ta converter melts (T b > 3270 K) 

around 2280 K. It is clear that effective cooling of the 
converter is essential at these high currents. 

B. Variation with converter thickness 

Figure [15] shows plots of (a) the maximum current I max 
and (b) the maximum yield as a function of converter 
thickness for a range of electron energies between 5 and 
100 MeV. The simulations suggest that if the electron 
current from the accelerator is not a limiting factor then 
increased positron production is possible if thicker than 
optimal converters are used at low energies (E e < 20 
MeV). At lower energy, increasing the thickness beyond 
optimal has little impact on the fractional energy ab- 
sorption (already near saturation), hence the maximum 
possible electron current increases and negates the slight 
reduction in j s . Conversely, for high electron energies 



E d I e = P w + P r = + 2ea( f (T(r) 4 - T^rdr) 

(0.5 + m(r 2 /ri)) J 

current I max is given by, 

_ C(T m - T w ) + D[Q (T(r) 4 - T 4 )2^dr] 

-'-max ' ^ 



TABLE II. Calculated values of the normalized slow positron 
production probability, 7jv, and energy deposited in a tanta- 
lum converter, Ed, for normal incidence electron irradiation 
on converters of optimal thickness, l opt . Values of the of the 
slow positron beam intensity, Y , and temperature, Tb, of the 
converter at an electron current of 10 /xA and 200 /iA are 
also included. Numbers in square brackets indicate order of 
magnitude. 



10 nk 200 fj,A 



E e 


7jv 


E d 


Y 


Tb 


Y 


T b 


(MeV) 


(e+/10 8 e-; 


(MeV/e~) 


(e+/s) 


(K) 


(e+/s) 


(K) 


5 


0.35 


2.44 


2.10[5] 


434 


4.19[6] 


2277 


10 


2.78 


5.55 


1.66 [6] 


389 


3.33 [7] 


2296 


15 


6.91 


7.83 


4.14[6] 


394 


8.28 [7] 


2234 


20 


12.5 


10.0 


7.47[6] 


399 


1.49 [8] 


2239 


40 


39.2 


16.6 


2. 35 [7] 


411 


4.70 [8] 


2297 


70 


83.5 


24.1 


5. 00 [7] 


423 


1.00 [9] 


2525 


100 


128.6 


31.5 


7.68[7] 


431 


1.50 [9] 


2790 



the optimum converter thickness for positron production, 
lopt, occurs at a point on the fractional energy absorp- 
tion curve (figure [TBla) - ) below the saturation level. This 
means that increasing the thickness beyond optimal val- 
ues also increases the energy deposition and the hence 
lowers the maximum possible electron current. 



VI. SUMMARY 

In conclusion, the present Monte Carlo simulations 
of the convertor and moderator provide an estimate of 
the slow positron yield as a function of electron energy 
and current provided by an electron accelerator in the 
range 5 < E e < 100 MeV. Although there is an expo- 
nential decrease in the positron production probability 
from the electron-positron converter as the electron en- 
ergy is reduced, the lower energy distribution of the emit- 
ted positrons results in slightly higher moderation effi- 
ciencies. Intense slow positron beams can be generated if 
high electron currents are used, for example at 10 MeV 
and 200 [iA an intensity of 3.3 x 10 7 e + /s is predicted, 
similar to the typical intensities achieved with existing 
high energy, low current facilities. 

The limiting factor for high current operation may not 
be the accelerator performance but rather the power de- 
posited in the converter. For low energy, high current 
accelerators, the maximum slow positron yield may be in- 
creased by using converters of greater than optimal thick- 
ness. The design of the converter and moderator assem- 
bly for any low energy, high current accelerator will need 
to consider the beam heating of both the converter and 
moderator, rotating targets will improve the maximum 
possible beam power but their design and implementa- 
tion would add further complexity to the system. 
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FIG. 15. (Color online) (a) Maximum electron current, Imax, 
based on equation [7] and (b) maximum positron yield, Y max , 
based on the product of the "f s and I ma x as a function of 
converter thickness, for electrons normally incident on a 
tantalum converter. 
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